Numerical evidence of the spin-chirality decoupling in the three-dimensional 

Heisenberg spin glass 



oo : 
o . 
o ■ 

(N ; 

s : 
<: 

in ■ 



I 



i 

c 

o 
o 



> 

00 
(N 

m 
co 

od 

o 

oo 
O 



Dao Xuan Viet and Hikaru Kawamura 
Department of Earth and Space Science, Faculty of Science, Osaka University, Toyonaka 560-0043, Japan 

(Dated: August 25, 2008) 

Ordering of the three-dimensional Heisenberg spin glass with Gaussian coupling is studied by 
extensive Monte Carlo simulations. The model undergoes successive chiral-glass and spin-glass 
transitions at nonzero temperatures Tea > Tsc > 0, exhibiting the spin-chirality decoupling. 

PACS numbers: 



The issue of the spin-glass (SG) ordering has been 
studied quite extensively for years, and continued to give 
an impact on surrounding areas. Meanwhile, the origi- 
nal problem of the magnetic ordering of typical spin-glass 
magnets e.g., canonical SG, still remains to be elusive 
As magnetic interactions in many SG materials are nearly 
isotropic, it is important to elucidate the ordering prop- 
erties of the three-dimensional (3D) isotropic Heisenberg 
SG. Although earlier numerical studies suggested that 
the 3D Heisenberg SG exhibited only a T = transition 
0,0], one of the present authors (H.K.) suggested that 
the model might exhibit a finite-temperature transition 
in its chiral sector 0. Chirality is a multispin variable 
representing the handedness of the noncollinear or non- 
coplanar structures induced by frustration. It has sub- 
sequently been suggested that, in the ordering of the 3D 
Heisenberg SG, the chirality was "decoupled" from the 
spin, the chiral-glass (CG) order taking place at a tem- 
perature higher than the SG order, T C g > T S g IEIE0I- 
Based on such a spin-chirality decoupling picture of the 
3D isotropic Heisenberg SG, a chirality scenario of exper- 
imental SG transition was proposed [J, : According to 
this scenario, the chirality is a hidden order parameter of 
SG transition. Real SG transition of weakly anisotropic 
SG magnets is then a "disguised" CG transition, where 
the chirality is mixed into the spin sector via a weak ran- 
dom magnetic anisotropy. 

Although consensus now seems to appear in recent nu- 
merical studies that the 3D HeisenbergSG indeed ex- 
hibits a finite-temperature transition 0,111, IE S B U, EH 
[Tl| . the nature of the transition, especially whether the 
model really exhibits the spin-chirality decoupling, is still 
under hot debate. Obviously, it is crucially important to 
clarify the ordering of the 3D Heisenberg SG model. 

The present situation, however, is not completely satis- 
factory. Mentioning some of the recent numerical works: 
By simulating the model of modest lattice sizes L < 20 {L 
being the linear dimension) but with rather small number 
of samples of N s — 32 (for their largest L), Hukushima 
and Kawamura presented support for the spin-chirality 
decoupling [f|. By contrast, Lee and Young claimed on 
the basis of their data of the correlation- length ratios £/L 
that the spin and the chirality order at a common tem- 
perature, thus no spin-chirality decoupling [jj [ll[ . Their 
data, however, suffers from either small lattice sizes of 



only L < 12 |9[ or small number of samples of N s = 56 
[ill ]. Campos et al simulated the same model to much 
larger lattices L = 32 with larger number of samples 
N s = 1,000, but no data below T g [Io| . Campos et al 
claimed that the chiral and spin sectors undergo simulta- 
neously a Kosterlitz-Thouless (KT) transition with mas- 
sive logarithmic corrections. This interpretation, how- 
ever, was criticized in Ref . [l2l] . 

Under such circumstances, we perform here a large- 
scale Monte Carlo simulation of the 3D Heisenberg SG 
in order to shed further light on the nature of its spin 
and chirality ordering. We exceed the previous simula- 
tions by simulating the system as large as L = 32 to 
temperatures considerably lower than T g for large num- 
ber of samples of order N s — 10 3 . Note that none of the 
previous simulations satisfied all these criteria simulta- 
neously. More importantly, we calculate several indepen- 
dent physical quantities including the correlation-length 
ratios, the Binder ratios and the glass order parameters, 
trying to draw consistent picture from these independent 
quantities, whereas Refs. [3, [13, [HI concentrated almost 
exclusively on the correlation-length ratio. Our simula- 
tion then enabled us to conclude that the SG transition 
occurs at a nonzero temperature which is located about 
15% below the CG transition temperature. Thus, the 
3D Heisenberg SG certainly exhibits the spin-chirality 
decoupling. 

The model is the isotropic classical Heisenberg model 
on a 3D simple cubic lattice with the nearest-neighbor 
Gaussian coupling. The Hamiltonian is given by 



« = - E J > 



i j Si ' Sj 



(i) 



where Si = (Sf,Sf,Sf) is a three-component unit vec- 
tor at the i-the site, and the < ij > sum is taken over 
nearest-neighbor pairs. The coupling Jy are random 
Gaussian variables with zero mean and the variance J 2 . 

The local chirality at the i-th site and in the /x-th direc- 
tion Xifj, m ay be defined for three neighboring Heisenberg 
spins by the scalar 



(2) 



where e M (// = x, y, z) denotes a unit vector along the 
/i-th axis. There are in total 3Af local chiral variables. 



2 



The lattice contains N = L 3 sites with L = 6, 8, 12, 
16, 24, 32 with periodic boundary conditions. Sample av- 
erage is taken over 2000 (L = 6, 8, 12), 1000 (L = 16, 24) 
and 800 (L = 32) bond realizations. To facilitate effi- 
cient thermalization, we employ the single-spin-flip heat- 
bath and over-relaxation method (Iol | , combined with the 
temperature-exchange technique. Over-relaxation sweep 
is repeated L times per every heat-bath sweep. 

Care is taken to make sure that the system is fully 
equilibrated. Equilibration is checked by the following 
procedures. First, we monitor the system to travel back 
and forth many times along the temperature axis during 
the temperature-exchange process (typically more than 
10 times) between the maximum and minimum temper- 
atures, while we also check that the relaxation due to 
the single-spin-flip updating is fast enough at the high- 
est temperature. This guarantees that different parts 
of the phase space are sampled in each "cycle" of the 
temperature-exchange run. Second, we follow Ref. (13} 
and check the equality expected to hold for the model 
with Gaussian couping. Third, we check the stability of 
the results against at least three times longer runs for a 
subset of samples. Fourth, we compare the data of the 
correlation-length ratios with the recent data by other 
authors in the temperature range where common data 
are available [13, El- Error bars are estimated by the 
sample-to-sample statistical fluctuations. 

We run two independent systems (1) and (2), and cal- 
culate a fc-dependent overlap. For the chirality, the k- 
dependent chiral overlap q x (k) is defined by the scalar, 



1 N 
= 3jy $J v <" 

i— 1 fi—x 7 y,z 



(1) (2) ik-f, 



(3) 



whereas, for the spin, it is defined by the tensor q a p(k) 
between the a and (3 components of the Heisenberg spin, 



N 



(1) c(, 2 )jk-fi 
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(a,/3 = x,y,z). (4) 



Note that both gcG and gsG are normalized so that they 
vanish in the high-temperature phase for L — > oo and 
gives unity in the nondegenerate ground state as expected 
for the present Gaussian coupling. 

The finite-size correlation lengths are given by Q 



l 



/ [(g(Q) 2 )] 

2sin(fc m /2)y [(g(fc m )2)] 



1, 



(9) 



for each case of the chirality and the spin, £cg and ^g, 
where k m — (2tt/L,0,0) with k m — \km\, and the \i- 
direction in Eq.(2) is taken here being parallel with k. 
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The CG and SG order parameters are defined by the 
second moment of the k = component of the overlap, 



* (fc) 2 = 51 \q a p(k) 

cx,/3—x,y,z 



(•5) 



(6) 



where (• • • } represents the thermal average and [• • • ] the 
average over the bond disorder. The chiral and spin 
Binder ratios are defined by 
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FIG. 1: The temperature and size dependence of the 
correlation- length ratio for the chirality (a), and for the spin 
(b). The arrow indicates the bulk chiral-glass transition point. 

In Fig.l, we show the correlation-length ratios for the 
chirality ^cg/L (a), and for the spin £sg/L (b). While 
the chiral £cg/L curves cross at temperatures which are 
only weakly L-dependent, the spin ^sg/L curves cross 
at progressively lower temperatures as L increases. The 
£/L data are compared with the data by other authors as 
follows: Our data for £/L are in full agreement with those 
of Ref.fToj within statistical error bars over the narrow 
and relatively high-temperature range covered by their 
data. The data of Ref . [Il|] for their largest L (on which 
their claim for "merging" was based) are lower than our 
present ones and those of Ref. [l(| by about 5 to 6 of our 
a units; this may be a purely statistical effect in view of 
the limited number of samples measured in Ref. (TlT| . 
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FIG. 2: The (inverse) size dependence of the crossing tem- 
peratures of £cg/L and £sg/L, the dip temperature Tdi p and 
the crossing temperature T cr oss of gca- From the linear ex- 
trapolation of the data, the sin-glass and chiral-glass transi- 
tion temperatures are estimated as Tcg = 0.145 ± 0.004 and 
Tsc = 0.120 ± 0.006. The inset exhibits a wider range. 



To estimate the bulk CG and SG transition tempera- 
tures quantitatively, we plot in Fig. 2 the crossing temper- 
ature of £c7G / L and ^sg I L for pairs of successive L values 
versus 1/L av , where L av is a mean of the two sizes. The 
data show an almost linear 1 / L at ,-dependence. The chiral 
crossing temperature exhibits a weaker size dependence, 
and is extrapolated to Tcg = 0.145±0.004 (in units of J), 
while the spin crossing temperature exhibits a stronger 
size dependence extrapolated to Tsg = 0.120 ± 0.006. 
Hence, Tsg is lower than Tcg by about 15%. For our 
data, we also tried a KT scaling with a logarithmic 
correction as was done in Ref . [l(| • KT-like scaling can be 
ruled out, however, when our lower temperature tag IT 
and £sg/L data, which were not available to Ref. [10f| . 
are included. This is particularly clear for ^cg/L data 
where the curves are not "merging" [l(| nor "marginal" 
[lH . but splay out below Tcg- 

The Binder ratios are shown in Fig. 3 for the chiral- 
ity (a), and for the spin (b). The chiral Binder ratio 
gcG exhibits a negative dip which deepens with increas- 
ing L. The data of different L cross on the negative 
side of gcc- These features indicate a finite-temperature 
transition in the chiral sector. By extrapolating the dip 
temperature Tdi P and the crossing temperature T cross for 
pairs of successive L values to L — oo, Tcg is estimated 
to be T CG = 0.145 ± 0.005: See Fig.2. The estimated 
Tcg is fully consistent with the one estimated above 
from {;cg/L. The peculiar form of gcc with a negative 
dip is consistent with the occurrence of a one-step-like 
replica-symmetry breaking as suggested by Hukushima 
and Kawamura Indeed, the calculated chiral- 

overlap distribution below Tcg ( n °t shown here) devel- 
ops a prominent central peak at q x = similar to the one 
reported in Refs. [!, 0] . 

By contrast, the corresponding spin Binder ratio gsc 
does not exhibit crossing nor merging in the temperature 
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FIG. 3: The temperature and size dependence of the Binder 
ratio for the chirality (a), and for the spin (b). The arrow 
indicates the bulk chiral-glass transition point. 

range studied, suggesting that the SG transition temper- 
ature, if any, occurs below T ~ 0.12. Meanwhile, as the 
size L is increased, gsc develops more and more singular 
form at low temperature, indicating that the associated 
overlap distribution significantly changes its shape at low 
temperature. If one recalls the fact that gsc takes a value 
unity at T = 0, gsG is expected to develop a negative dip 
at a lower T (of IS 0.12) accompanied by an upturn to- 
ward T — 0. This feature of gsc suggests the occurrence 
of a SG transition at a nonzero temperature, Tsg ~ 0.12. 

Pixley and Young recently criticized that the Binder 
ratio might not be an appropriate quantity in studying 
the ordering of vector SG, arguing that the large num- 
ber of the order-parameter components (3 2 = 9 in the 
Heisenberg SG) might lead to a trivial Gaussian distri- 
bution even below T g [3]. In fact, mere the large number 
of order-parameter components does not lead to such a 
trivial behavior: For example, in a simple 3D 0(n) ferro- 
magnet with n — 6 fl~5| and 10 (l6| . the spin Binder ratio 
exhibits at T c a clear crossing behavior characteristic of 
the standard long-range ordered phase, quite different 
from the one of Fig. 3(b), presenting counter-examples to 
the criticism of Ref. fl4|] . Hence, the peculiar behavior of 
gsc observed here in Fig. 3(b) should be regarded as a 
manifestation of essential features of the SG ordering. 

In Fig. 4, we show the size dependence of the CG and 
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FIG. 4: Size dependence of the chiral-glass order parameter 
Iqq (a), and the spin-glass order parameter q)g<L (b). Straight 
lines in the figures are drawn by fitting the three data points 
of smaller sizes, L — 6, 8 and 12. 



SG order parameters on a log- log plot. As can be seen 

(2) 

from Fig. (a), q CG exhibits a linear behavior at a temper- 
ature T = 0.148, a downward curvature characteristic of 
a disordered phase at higher T, and an upward curva- 
ture characteristic of a long-rage ordered state at lower 

(2) 

T. Thus, the data of q C Q are consistent with a CG tran- 
sition occurring at Tec = 0.148±0.005, consistently with 
the results of £cg/L and gcc- 



(2) 

The SG order parameter q S( j exhibits a significantly 
different behavior, i.e., it exhibits a downward curvature 
characteristic of a disordered state at T = 0.148 ~ Tea, 
or even at T = 0.133 < Tcg- At our lowest temperature 
T/J = 0.121 where we could equilibrate only smaller lat- 
tices of L < 16, the data exhibit a near linear behavior up 
to L = 16, although it is not clear whether this linear be- 

(2) 

havior extends to larger L. Thus, our data of qsc(L) are 
consistent with a SG transition occurring at T$g ~ 0.13, 
whereas, from the present data of q s 2 Q only, we cannot 
rule out the possibility that T$g is significantly lower 
than this. 

All the physical quantities studied here, including the 
correlation-length ratio, the Binder ratio and the glass or- 
der parameter, consistently indicate that the 3D Heisen- 
berg SG with Gaussian coupling exhibits successive CG 
and SG transitions at T C g = 0.145 ± 0.004 and at 
Tsg ~ 0.12. The SG order sets in at a temperature 
at least about 15% below the CG order, hence, the oc- 
currence of the spin-chirality decoupling. One may feel 
that the relative distance between Tcg and Tsg is not so 
large, but, in fact, it is a sizable difference, much larger 
than the one observed in other systems exhibiting the 
spin-chirality decoupling, e.g., the 2D regular frustrated 
AT model where the difference is known to be less than 
1% [H El. While the SG order in the 3D Heisenberg 
SG occurs at nonzero temperature, as is consistent with 
the recent numerical works 

H & M EH, it should be 
stressed that whether Tsg is zero or nonzero is irrele- 
vant to the chirality scenario of Refs.fj, |7| as long as the 
spin-chirality decoupling occurs, i.e., Tsc < Tec- Our 
present result then supports the chirality scenario of SG 
transition. 
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